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20.  Abstract  (Contd) 

>  characteristics  of  the  scattering  from  the  surface  depend  on  the  degree  of 
roughness..  Tee  first  is  the  changes  in  form  of  the  mathematical  represen¬ 
tation  forjji0  as  the  degree  of  roughness  is  varied.  The  next  topic  is  con¬ 
cerned  with  alternative  definitions  of  the  glistening  surface  (the  part  of  the 
rough  surface  from  which  diffuse  power  is  scattered  into  the  receiver).  It 
is  shown  that,  for  some  parameter  ranges,  the  usual  definition  of  glistening 
surface  as  given  by  Beckmann  and  Spizzichino  is  not  adequate.  In  the  final 
topic,  the  effect  of  the  moisture  content  of  the  soil  on  surface  scattering  as 
a  function  of  roughness  is  discussed.  A  wide  range  of  results  is  presented 
in  each  of  these  areas.  The  implications  of  these  findings  for  the  prediction 
of  system  performance  are  discussed.  _ 
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Further  Considerations  in  Models 
of  Rough  Surface  Scattering 


1.  INTRODUCTION 

The  scattering  of  electromagnetic  signals  from  rough  terrain  depends  on  the 

characteristics  of  the  surface.  The  scattering  surface  can  be  defined  in  terms  of 

the  statistical  distribution  of  the  heights,  their  degree  of  correlation,  the  variance 
2 

of  the  heights,  a  ,  and  the  complex  dielectric  constant  appropriate  to  the  type  of 

1  -5 

terrain.  Theoretical  models  describing  rough  surface  scattering  relate  these 
parameters  to  the  normalized  cross  section  of  the  terrain.  This  report  is  con¬ 
cerned  with  several  factors  that  can  affect  the  description  of  these  relationships. 


(Received  for  publication  7  January  1083) 
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At  RADI'  EEC  there  is  a  continuing  interest  in  impri'V':  modeling  of  tin  -  scat¬ 
tering  from  rough  surfaces.  The  investigations  'oncentrate  on  two  main  themes: 
the  description  of  the  surface  and  the  electromagnetic  interaction.  Thr  range  of 
models  can  be  seen  by  the  descriptions  in  earlier  reporting.  *'  l(  The  details  of 
the  models  will  be  presented  in  later  sections  and  in  the  appendices.  It  should 
just  be  noted  here  that  tne  models  are  comprehensive  and  have  the  capability  of 
calculating  the  coherent  (specular  plus  direct)  and  incoherent  multipath  power 
reaching  a  monopulse  receiver  from  a  beacon  located  over  rough  terrain.  Dis¬ 
cussions  of  coherence  and  the  methods. for  calculating  the  coherent  and  incoherent 
power  are  given  in  the  appendices.  One  additional  point  is  that  in  these  models  we 
use  the  concept  of  "glistening  surface"  (for  a  given  set  of  antenna  positions,  that 
part  of  the  rough  surface  that  reflects  a  significant  amount  of  the  transmitted 
ene  rgy  into  the  receiver)  as  described  by  Beckmann  and  Spizzichino.  1  A  geometri¬ 
cal  representation  of  the  scattering  from  the  surface  is  shown  in  figure  1. 


«ECT  S1ON.V 


'  'L.V  V  -iPAm 


:vc.r?f  PECt  .£« 


SIFF.SE  V  ,\T  .'V. 


X  S'CNiNd  'V  PEACE  -U  The  Pl  ANE  OF  Th:  EARTH 

Figure  1.  Geometrical  Representation  of  the  Scattering  From  Rough 
Surfaces 


3.  Papa,  R  .J. ,  and  Lennon,  J.  F.  (1980)  Electromagnetic  scattering  from  rough 
surfaces  based  on  statistical  characterization  of  the  terrain.  International 
Radio  Science  Symposium,  (URSI),  June  1980),  Quebec,  Canada. 

1.  Papa,  R.J.,  Lennon,  J.F.,  and  Taylor,  R.  L.  (1980)  Electromagnetic  Wave 
Scattering  From  Rough  Terrain,  RADC-TR -80-300,  AD-A-n>i839. 


In  this  report  we  will  be  addressing  three  different  topics  related  to  the  scat¬ 
tering  from  rough  surfaces  and  will  show  how  all  of  them  are  sensitive  to  the 
actual  degree  of  roughness.  The  first  is  the  determination  of  the  form  for  the 
normalized  cross  section  of  the  surface,  a0.  The  usual  form  is  based  on  physical 
optics  and  the  asymptotic  evaluation  of  integrals  representing  a-  in  the  limit  of 
very  rough  surfaces  for  bivariate  Gaussian  and  exponential  surface-height  distribu¬ 
tions.  In  the  exponential  case  the  uncorrelated  variates  are  not  independent  and, 
for  intermediate  levels  of  roughness,  this  leads  to  a  form  for  a®  different  from  the 
usual  result.  The  next  topic  is  that  of  the  glistening  surface,  ^ar  some  parameter 
ranges,  the  classical  definition  of  the  glistening  surface  as  given  by  Beckmann  and 
Spizzichino  is  not  adequate.  Their  definition  does  not  account  for  the  angular 
dependence  of  the  scattering  matrix  in  the  expression  for  a®  (the  normalized  cross 
section  of  the  rough  surface).  In  addition,  it  may  suppress  contributions  from 
areas  close  to  the  antennas.  The  final  roughness  dependent  effect  discussed  here 
is  how  the  scattering  from  a  given  type  of  surface  is  effected  by  the  amount  of 
moisture  present  in  the  soil.  To  represent  this,  the  dielectric  constant  of  the  sur¬ 
face,  for  a  range  of  conditions,  is  chosen  to  correspond  to  arid,  normal,  and  very 
wet  soil.  The  effect  is  again  sensitive  to  the  degree  of  roughness. 

The  results  of  this  study  can  have  implications  for  the  estimation  of  radar  sys¬ 
tem  performance.  Two  aspects  arc  important,  first,  there  is  the  question  of 
establishing  system  design  criteria  that  are  cost  effective.  Second,  the  assess¬ 
ment  of  system  performance  variability  is  related  to  the  discussions  of  the  effect 
of  surface  roughness  an  moisture  dependent  scattering. 

2.  Till.  NOK MAI.IZKI)  CROSS  SUCTION  01  RANDOM  ROUGH  SliRFACKS 

In  this  section  we  will  discuss  the  basic  theoi  etical  considerations  involved  in 
the  determination  of  tne  rough  surface  scattering  cross-section.  Particular  con¬ 
sideration  will  be  given  to  limitations  of  the  theories  and  the  effects  of  the  statis¬ 
tical  properties  assumed  to  describe  the  surface  h°ight  distribution  of  the  scatter¬ 
ing  area. 

2.1  >'iii"lf’-Si'alt’-ol-K<>i!"lim,*.«  Theories 

Two  major  theoretical  formulations  for  the  normalized  cross  section,  the 
physical  optics  model  and  the  perturbation  theory  model,  will  be  considered 
in i  hiding  the  delineation  of  conditions  under  which  each  applies. 


•J.  1.  1  PHYSICAL  OPTICS 


JO 

In  thr  physical  optics  approximation,  ' "  the  radius  of  curvature  of  the  irregu¬ 
larities  is  assumed  to  be  large  compared  to  a  wavelength  (II  _  »  \).  The  recent 

4  ° 

work  by  Brown  has  shown  that  this  condition  is  not,  in  itsel'1,  sufficient  for  phys¬ 
ical  optics  to  be  applicable.  An  additional  cjnstraint  is  required,  such  as  the 
surface  slopes  must  also  be  small  (T  »  o),  where  T  =  surface  correlation  length 
and  o  -  standard  deviation  in  surface  height.  Also,  the  total  field  on  the  rough 
surface  may  be  expressed  as  the  sum  of  the  incident  field  plus  the  scattered  field. 
Because  the  radii  of  curvature  are  assumed  to  be  large  compared  to  a  wavelength, 
the  scattered  field  is  expressed  as  the  incident  field  multiplied  by  the  Fresnel 
plane  wave  reflection  coefficient.  Then,  the  normalized  scattering  cross  section 
(a0)  for  the  surface  is  derived  by  substituting  this  approximation  for  the  total  field 
on  the  surface  into  the  Kirchhoff  integral  (scalar  or  vector  form)  expression  for 
the  scattered  field.  The  expression  that  is  derived  for  the  normalized  cross  sec¬ 
tion  o’  of  the  rough  surface  depends  upon  the  probability  distribution  function 
(PDF)  of  the  urface  heights,  the  surface  correlation  function,  standard  deviation 
in  surface  height,  surface  correlation  length,  and  complex  dielectric  constant  of 
the  surface. 

2.1.2  PERTURBATION  THEORY 

The  second  conventional  method  for  deriving  an  expression  for  a*  with  a  single 
scale  of  surface  roughness  is  the  Rayleigh-Rice  or  perturbation  method.  Here,  it 
is  assumed  that  the  scattered  or  transmitted  fields  from  or  across  a  rough  surface 
can  be  represented  by  a  superposition  of  plane  waves  propagating  away  from  the 
boundary  with  unknown  amplitudes  (Rice,  Peake,  and  Ruck  et  al  ).  Also,  it  is 
assumed  that  the  standard  deviation  in  surface  height  is  small  compared  to  a  wave¬ 
length  and  that  the  surface  slopes  are  small.  Then,  the  scattered  fields  are  ex¬ 
panded  in  a  perturbation  series.  The  unknown  scattered-field  amplitudes  are 
determined  by  requiring  each  order  of  perturbation  to  satisfy  Maxwell’s  equations 
and  the  boundary  conditions.  The  expression  that  is  derived  for  a*  is  proportional 
to  the  surface  height  spectral  density,  the  fourth  power  of  the  wavenumber,  and  the 
fourth  power  of  the  sui  face -height  standard  deviation. 

2.2  Multiple  Scale*  of  Roughnrs* 

There  have  been  several  models  developed  to  describe  wave  scattering  from 
a  surface  with  two  or  more  scales  of  roughness.  Some  two-scale  models  have  been 


9.  Rice,  S.D.  (1951)  Reflection  of  electromagnetic  waves  by  slightly  rough 

surfaces.  The  Theory  of  Electromagnetic  Waves,  M.  Kline,  Ed.,  Wiley, 
New  York.  . 
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derived  by  fcssuunng  that  tr.e  small  scale  of  roughness  obeys  the  criteria  for  the 
perturbation  solution  that  is  superimposed  on  a  lai  ge  scale  of  roughness  tnat  obeys 
tlie  physical  optics  ci  tterion  (Wright*l)  and  Fuks**).  Brown*”  has  obtained  an 
expression  for  e*  with  two  scales  of  roughness  using  a  perturbation  technique  first 
derived  by  Burroughs.  Bahar  nas  used  the  full  wave  approach  to  derive  ex¬ 
pressions  for  o*  to  describe  the  scattering  of  em  waves  from  a  rough  surface  with 
more  than  one  scale  of  roughness.  The  solutions  have  been  shown  to  reduce  to 
the  physical  optics  result  or  the  perturbation  result  in  appropriate  limits. 

2.1  Surface  Height  Statiatii 

2.3.  1  GAUSSIAN  HEIGHT  VARIATES 

Most  authors  have  considered  that  the  random  rough  surface  from  which  the 
scattering  occurs  is  Gaussian  in  nature.  There  are  two  main  reasons  for  this 
assumption.  In  the  first  case,  there  are  a  number  of  instances  in  which  agree¬ 
ment  with  experiment  has  been  obtained  with  this  assumption.  Secondly,  this 
assumption  is  attractive  because  of  the  statistical  properties  of  Gaussian  distribu¬ 
tions.  For  instance  ,  sums  of  Gaussian  variates  are  also  Gaussian.  Another 
example  of  perhaps  more  interest  is  the  property  that  once  the  variates  become 
uncorrelated,  their  respective  marginal  densities  are  independent.  Both  of  these 
aspects  are  important  in  the  mathematical  analyses  for  the  scattering  cross  sec¬ 
tion  of  the  surface. 

2.3.2  NON -GAUSSIAN  INDEPENDENT  VARIATES 

The  first  serious  effort  to  investigate  em  -wave  scattering  from  non-Gaussian 

15 

rough  surfaces  is  presented  in  the  report  by  Barrick.  Here,  expressions  for  a* 
are  derived  under  the  assumptions  of  physical  optics  for  a  modified  Bessel  joint 
PDF  and  a  Gaussian  joint  PDF.  Both  a  Gaussian  and  an  exponential  function  are 

10.  Wright,  J.W.  (]9G(i)  Baekscattering  from  capillary  waves  with  application  to 

sea  clutter,  IEEE  Trans.  Antennas  Propag.  AP-14:749. 

11.  Fuks,  I.  M.  (19GG)  Contribution  to  the  theory  of  radio  wave  scattering  on  the 

perturbed  sea  surface,  Izv«’stia  Vyshikh  Uchebnykh  Zavedenii  Kadiufizika 
3:87(5.  ~  *  . "  "  ~  ~ 

12.  Brown,  G.S.  (I:>73)  Baekscattering  from  a  Gaussian  distributed  perfectly 

conducting  rough  surface,  IEEE  Trans.  Antennas  Propag.  AP-2K:47J. 

13.  Burrows,  M.  L.  (1937)  A  reformulated  boundary  perturbation  theory  in  elec¬ 

tromagnetism  and  its  application  to  a  sphere.  Can.  J,  Phys.  43: 1729. 

14.  Bahar,  E.  ( 1  t )  -Scattering  cross  sections  for  random  rough  surfaces:  Full 

wave  analysis,  Radio  Sci.  PB331-341. 

1 5 .  Barrick,  D.  1-  .  (jnc~)  A  More  Exact  Theory  of  Baekscattering  From  Statis¬ 

tically  Rou„h  Surfaces,  T>hio  State  University  Research  Foundation, 

Report  IT! 3 -13. 
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.issumt  !  I'm  the  :~ui  f:u  >  -in  iglit  c.  r  relation  I'm  -’ion.  ii.-t  r  •  i-  is....  use  :  ,i  l.-i  i.i  i  n 
the  modified  Bessel  joint  FIJI-  that  was  first  proposed  by  4 >tt.  1  If  Z  and  rep- 
ri'si'nt  surface  heights  at  two  -list met  points  on- the  surface,  ’hi  ti  tin-  modified 
Bessel-joint  PDF  is  given  by 


PB’V  '•o  ^ 


Vil  z  -  cz..r 

"U3  i  Z„  | 

)  i42!  i/.,  -  eZ,i  K, 

K1 

IT 

when  <i  is  tin  standard  deviation  in  surfaei  height,  e  is  tin  . •om-tatmn  function, 
and  K,  ;s  the  modified  Bessel  fun-  tion  of  the  first  kind  of  older  one.  It  should  be 

noted  that  -.'.hen  t  =  (x^  -  x., )“  +  (y  -  y.,)“  —  *  and  <  (t>  -  then 

Pj^(/.  j ,  )  -  [>  (/.  )p.,(Z.,).  This  means  that  decorrelation  implies  statistical  inde¬ 

pendence  for  this  PDF.  However,  the  two  variates  are  not  described  by  identh  al 
marginal  probability  densities  in  this  ease  and  tin-  meaning  of  tins  result  is  some¬ 
what  unclear. 

In  1P73,  Beckmann*1  investigated  t'ra- "ttve  scattering  from  rough  non -( biussian 
surfaces  using  the  assumptions  of  physical  optics.  The  method  assumes  that  the 
joint  PDF  pfZ^.Z^.c)  can  be  .erived  from  the  marginal  density  pfz)  and  known 
correlation  function  e(f).  The  technique  uses  the  properties  of  classical  sets  of 
polynomials  (llermite,  generalized  Laguerre,  and  Jacobi),  which  are  orthogonal 
over  some  interval  with  respect  to  some  weighting  function.  Another  assumption 
is  that  the  joint  PDF  p(Z^,  Z^,  c)  -  as  c(r  —  a)  —  o,  that  is,  deror- 

t  elation  leads  to  statistical  independence. 

2.3.3  NON -GAUSSIAN  DFPKNDFNT  VAKIATFS 

The  transformation  of  joint  surface-height  probability  densities  into  sets  of 
uncorrelated  variates  does  not  always  lead  to  statistical  independence.  In  Buck 

9 

et  al,“  Barrtek  introduced  the  bivariate  exponential  surface  height  PDF; 


pe(5VV 


2 ,  9,772 

2m’  (1  -  cl1'" 


)  exp 

[•( 

Z“  -  2c^j/'2  +  ^2 
(l/3)eT(l  -  c2) 


1/2 


The  uncorrelated  result  for  this  case  is  not  a  pair  of  independent  variates,  and 
this  leads  to  complications  in  the  theoretical  development  of  o®.  The  authors  have 


lb.  Ott,  K.H.  (1PP7)  A  Theoretical  Model  for  Scattering  Prom  Hough  Surfaces 
With  Applications  to  the  Moon  and  Sen,  Ohio  State  University  Research 
Foundation,  Report  - 1388-1. 

17.  Beckmann,  P.  (lf’73)  Scattering  by  non-Gaussian  surfaces,  IF  UK  Trans. 
Antennas  Fropag.  AP-2  1 :  p;f». 


fl 
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used  variations  of  this  distribution  in  several  previous  studies.  ’  *  As  a  re¬ 

sult  of  applying  hypothesis  testing  procedures  (see  Appendix  A)  it  was  found  that 
for  the  radar  site  {the  Discrete  Address  Beacon  System  in  Massachusetts 
(DABS)  }  used  for  experimental  comparison,  the  terrain  heights  fii  better  to  a 
bivariate  exponential  distribution  than  a  bivariate  Gaussian  distribution. 

Recently,  Brown"  has  considered  scattering  from  perfectly-conducting  ran¬ 
dom  surfaces  for  which  decorrelation  does  not  imply  statistical  independence.  By 
using  an  exact  theory  for  the  current  induced  on  the  surface  and  the  stochastic 
Fourier  transform  approach,  it  is  shown  that  the  incoherent  scattered  power  con¬ 
sists  of  two  parts.  The  first  part  corresponds  to  what  is  normally  termed  the  dif¬ 
fuse  power,  that  is,  the  incoherent  power  scattered  in  all  directions.  The  second 
part  is  incoherent  power  scattered  specifically  in  the  specular  direction. 


2.4  The  Condition*  For  Phy*it-al  Optica 

In  our  opinion,  the  physical  optics  models  of  rough  surface  scattering  are  the 
most  rigorous  of  all  models  to  date.  The  analyses  developed  herein  are  based  on 
the  limits  of  physical  optics  assumptions.  For  the  rough  surface  scattering  jdase 
we  can  examine  the  parameters  in  relation  to  those  conditions  for  which  physical 
optics  applies  In  this  analysis  the  rough  surface  heights  are  assumed  to  be  re¬ 
lated  by  Gaussian  correlation  functions. 

The  rough  surface  is  given  by  the  equation  z  =  ?{x,y).  For  an  isotropic  sur¬ 
face.  the  slopes  and  are  equal.  The  average  slope  is  given  by 


where  (  )  denotes  expectation  value.  F 

•>  ,  9 
2  -t-/T* 

relation  functional  is  R  *  o  e 

F  2  ')  1 

points,  t  *  ^(xj  -  x.,)  *  (y:  -  y.,)  J 


irst,  we  assume  that  the  surface  autoicor- 

j 

Here  t  is  again  the  separation  between 
The  square  of  the  radius  of  curvature. 


Rc  of  a  curve  is  given  by 


18.  Lennon,  J.F.,  and  Papa,  R.J.  (1980)  Statistical  Characterization  of  Rough 

Terrain.  RADC-TR-80-9,  AD  A087737T 

19.  Papa,  R.J.,  Lennon,  J.  F. ,  and  Taylor,  R.L.  (iP80)  Prediction  of  Electro¬ 

magnetic  Scattering  for  Rough  Terrain  Using  Statistical  Parameters  Derived 
From  Digitized  Topographic  Maps,  RADC-TR-80-28P,  ADA094104. 

20.  McGarty,  T.  P.  (1975)  The  Statistical  Characteristics  of  Diffuse  Multipath 

and  Its  Kffert  on  Antenna  Performance,  KSD-TK-75-145,  AD  A 009809. 

2  1.  Brown,  G.S.  (1982)  Scattering  from  a  class  of  randomly  rough  surfaces. 

Radio  Sci.  (to  be  published). 


The  expectation  value  of  It”  is  given  by 


>'A  =  /n  +  <g')~)3 

\  (z ")” 


To  simplify  tin*  evaluation  of  this  quantity  we  will  approximate  \R~  )  by  the  follow 

\  c/ 

imv  vxpn-ssion: 


/w 2  ‘ i±?l  '  =  U 


\(z  >7 


I’hr  average  of  (z')_  is  given  bv 


.,2  .  ! /as  _  3~tt (r  =  o) 


jV,x  » VI 


I  3T“ 


and  the  average  of  (z")“  is  given  by 


,  „,2  94H(t  =  0)  .  12a 

(z  )  = - 3 - - 7- 

I  3T  I  T 


Then,  an  approximate  expression  for  the  average  radius  of  curvature  is  given  by 


"c  *  -=-  [>  *  2«2/t2]‘ 
v  12  a 


We  have  already  discussed  the  fact  that  one  set  of  sufficient  conditions  for 

physical  optics  to  apply  is  R  »  A.  and  1  »  a/T.  If  we  consider  the  equations  for 

•>  ,  c 

It  and  a  new  condition  T”;.t  »  A  then,  these  relations,  along  with  the  relation 
T  o,  lend  to  tile  simpler  (in  the  sense  of  being  more  directly  seen)  set  of  suf- 
ficient  conditions  for  physical  opti<  s,  <t/T  <<  1  and  T  /a  A,  for  a  surface  with  a 
(■nussian  autocorrelation  function.  One  interesting  point  to  note  is  that  these  con¬ 
ditions  holt!  for  ail  values  nt  the  Rayleigh  roughness  parameter,  11  where 


L  =  1  ^  )  'a (cos  ♦  cos  W  ).  Here  A  =  em -wavelength  wavelength,  6^  =  angle  of 
incidence,  and  angle  of  scattering  measured  with  respect  to  surface  normal. 

2.5  o®  IMrriuiiulion 

In  this  chapter,  comparisons  will  be  made  between  the  exact  integral  represen¬ 
tation  for  a’  and  the  asymptotic  representation  for  o»  for  both  exponential  and 
Gaussian  POl-'s  as  a  function  of  tile  Rayleigh  parameter  E.  It  will  be  shown  that 
for  exponentially-distributed  surface  heights,  great  care  must  be  exercised  in  the 
use  of  the  expression  for  o®  as  given  in  the  book  by  Ruck  et  al.  In  Appendix  B, 
it  is  shown  how  expressions  for  cr®  may  be  derived  for  bivariate  exponential  PDFs 
under  the  assumptions  of  physical  optics.  For  large  Rayleigh  parameter  2,  a 
steepest  descent  evaluation  is  made  to  obtain  the  expression  for  a®  first  given  by 
Barrick  in  Ruck  et  al. 2  It  will  be  shown  in  the  results  of  this  cliapter  that  the 
asymptotic  expression  given  for  a®  is  not  accurate  unless  £  is  very  large.  This 
difficulty  arises  because  the  bivariate  exponential  PDF  does  not  reduce  to  inde¬ 
pendent  variates  Z^  and  Z0  as  t  =  £  (x^  -  x9)2  +  (v^  -  y2>"  j1^2  —  «,  and,  as  E 
increases,  the  normalized  cross  section  ct®  does  not  approach  zero  as  rapidly  as 
in  the  Gaussian  PDF  case. 

For  the  purpose  of  this  study,  the  behavior  of  a*  is  discussed  in  terms  of  the 
behavior  of  the  slope-dependent  statistical  quantity  J.  This  is  explained  in  detail 
in  Appendix  B.  The  formulations  are  developed  for  various  ranges  of  £.  The 
restrictions  and  regions  of  agreement  are  assessed  in  the  following  sections. 

2.  3.  1  DESCRIPTION  OF  TABULAR  RESULTS 

Tables  1  through  18  present  values  of  the  quantity  J  as  a  function  of  the  azi¬ 
muthal  scattering  angle,  PHISC,  for  both  Gaussian-distributed  surface  heights  and 
exponentially-distributed  surface  heights.  This  form  of  presentation  permits  us 
to  examine  whether  azimuthal  displacement  affects  agreement.  In  these  tables, 

£  is  the  Kavleigh  roughness  parameter  and  T  is  the  correlation  length  in  meters. 

2 

The  values  assumed  by  the  Rayleigh  parameter  are  given  as  £  =  95.864, 

2  0.  13,  10.  i),  1.  0,  and  0.  l  (for  exponential  surfaces).  The  correlation  length  is 

3  tn  and  15  m. 

For  Gaussian-distributed  surfaces,  SUM  refers  to  the  summation  representa¬ 
tion  for  J;  APPROX  refers  to  the  steepest  descent  evaluation  of  the  integral  repre¬ 
sentation  for  J  (valid  for  L  »  1);  the  INTEG  refers  to  the  integral  representation 
l"i  I.  For  the  exponentially-distributed  surfaces,  APPROX  refers  to  the  steepest 
descent  evaluation  of  the  integral  representation  for  J  (valid  for  £  »  l);  SING  INTEG 
refers  to  the  integral  representation  for  J;  and  POWER  EXPAN  refers  to  the  eval¬ 
uation  of  the  integral  representation  for  J  obtained  by  expanding  the  denominator  in 
a  power  series  (valid  for  £  ;  1).  These  are  all  derived  in  Appendix  B. 

1r' 


UtblP  1.  Values  of  J -integral  ns  a  Function  of  Azimuthal 
Scattering  Angle:  =  95.  804.  T  =  5.0,  Gaussian  Surface 


PH  ISC 


2. 005  • 

10‘14 

2.07  • 

10'14 

1.42  • 

10*14 

8.31  • 

10'1-’ 

4.  19  • 

10*1-’ 

1.31  • 

10" 1  ’ 

•  ».  7n 

to’1'5 

APPROX 

14  1009.34 

14 

951.70 

4  803.00 

- 

752.  47 
030.  9 1 
'  508.  70 

394.42 


I  NT  KG 

1019. 60 
900.30 
809.  02 
755.  77 
031.  78 
507.  77 
392.  48 


Table  2.  Values  of  J-integral  as  a  Function  of  Azi¬ 
muthal  Scattering  Angle:  £2  --  20.13,  T  =  5.0,  Gaus¬ 
sian  Surface 


PH  ISC 

SUM 

APPROX 

I NT  KG 

0.  5 

976.  13 

1009.34 

1002. 85 

1.0 

912.20 

951. 70 

995.29 

1.  5 

815.34 

803.  CO 

892.  72 

2.0 

697. 65 

752.  47 

707.  70 

2.  5 

572. 10 

630.  91 

033.  74 

3.0 

450. 19 

508.  70 

502. 92 

3.  5 

340. 41 

394.  42 

384. 26 

Table  3.  Values  of  J-integral  as  a  Function  of  Azimuthal 
Scattering  Angle:  £2  =  10.0,  T  =  5.0,  Gaussian  Surface 


SUM 

APPROX 

1132. 80 

1009.34 

1048. 11 

951.  70 

924.33 

803. 00 

780.01 

752.47 

032.  00 

030. 91 

493.38 

508. 70 

372. 81 

394.  42 

PH  ISC 


0.  5 
1.0 
1.  5 
2.0 
9 


I  NT  KG 


1131.37 
1049.21 
924.  153 
779.72 


493. 41 
371. 93 


Table  4.  Values  of  J-integral  as  a  Function  of  Azimuthal 
Scattering  Angle:  Z-  -  1.0,  T  =  5.0,  Gaussian  Surface 


PI  use 

SUM 

APPROX 

I NT  KG 

0.  5 

490. 02 

1009. 34 

0572. 220 

1.0 

400.  20 

951.  7  0 

-950,  11 

1.  5 

428.  49 

803. 00 

077.  40 

2.  0 

381.  00 

752. 47 

582  13 

2.5 

328.25 

030.  91 

-74.32 

3.  0 

271.28 

508.  70 

705,  20 

3.  5 

222.  70 

394.  42 

-290,  00 

Table  Values  of  J -integral  as  a  Function  of  A  z\ - 
muthr.l  Scattering  Angle:  £-  -  10.0,  T  =  15.0,  Gaus 
sian  Surface 


SIM 

APPROX 

INTEG 

8290.  55 

7752. 52 

8298.  12 

4430.  46 

4559. 19 

4431.  60 

1731.  12 

1893. 12 

1730.  97 

523.01 

551.  43 

527.  73 

131.54 

112.  93 

131.90 

27.74 

15.25 

27.77 

5.05 

1.  65 

4.  90 

PH  ISC 


Table  8.  Values  oi  J-intcgral  as  a  Function  of  Azi 
muthnl  Scattering  Angle:  =  1.0,  T  =  15.0,  Gaus 

sian  Surface 


PH  ISC 

SIM 

APPROX 

INTEG 

0.  5 

3850.28  , 

7752. 62 

9931.  88 

1.0 

2464.  43 

4569.  19 

1048.  07 

1.  5 

1221.  87 

1893. 12 

1470.  77 

2.  0 

506.  02 

551.  43 

707.39 

2.  5 

193. 80 

112. 93 

-208.  77 

3.0 

73.  71 

16.26 

564.  70 

3.  5 

27.  95 

1.  65 

-491.  41 

T;tblf  9.  Values  of  J-intt  steal  as  a  Function  of 
A zimuthn  1  Scattering  .Ant*!*':  11  ~  ••  95.110-1,  T  -  f>.  0 
Fxponent  in!  Surface 


Table  10.  Values  of  J-integral  as  a  Function  of 
Azimuthal  Scattering  Angle:  L*  =  20.  13,  T  =  5.0, 
[exponential  Surface 


Table  11.  Values  of  J-integral  as  a  Function  of 
Azimuthal  Scattering  Angle:  £2  =  10.0,  T  !  5,0, 
Exponential  Surface 


PH  ISC 

APPROX 

SING  INTEG 

0.5 

2188 

2395 

1.0 

1555 

1522 

1.  5 

1104 

1084 

2.0 

784 

769 

2.  5 

556 

526 

3.0 

395 

384 

3.  5 

280 

269 

Table  12.  Values  of  J-integral  as  a  Function  of  Azimuthal  Scatter¬ 
ing  Angle:  £2  =  1.0,  T  =  5.0,  Exponential  Surface 


PHISC 

APPROX 

SING  INTEG 

POWER  EXPAN 

0.  5 

2188 

10/8 

1009 

1.0 

1555 

3C8 

952 

1.  5 

1104 

398 

863 

2.  0 

784 

320 

752 

2.  5 

556 

273 

631 

3.0 

395 

234 

509 

L  3‘5 

280 

182 

394 

Table  13.  Values  of  J-intogral  as  a  Function  of  Azimuthal  Scatter¬ 
ing  Angle:  =  0. 1,  T  =  5.0,  Exponential  Surface 


PHISC 

APPROX 

SING  INTEG 

POWER  EX  PAN 

0.  5 

2188 

14 

1 

10 

1.0 

1555 

13 

10 

1.5 

1103 

11 

10 

2.0 

784 

10 

10 

2.5 

556 

9 

10 

3.0 

395 

9 

10 

3.  5 

280 

8 

9 

Table  14.  Values  of  J-integral  as  a  Function  of 
Azimuthal  Scattering  Angle:  £2  =  95.  864,  T  =  15.  0, 
Exponential  Surface 


PHISC 


APPROX 


SING  INTEG 


0.5 

1.0 


9881 

3544 

1269 

454 

162 

58 


9850 

3525 

1264 

454 

163 

58 


21 


21 


Table  15.  Values  of  J-integral  as  a  Function  of 
Azimuthal  Scattering  Angle:  E^  =  20.  13.  T  =  15.  0, 
Exponential  Surface 


APPROX 

SING  1NTEG 

9881 

9676 

3  544 

3450 

1269 

1247 

454 

453 

152 

155 

58 

60 

21 

22 

PH  ISC 


Table  16.  Values  of  J-integral  as  a  Function  of 
Azimuthal  Scattering  Angle:  E2  =  10.0,  T  =  15.  C 
Exponential  Surface 


PHISC 

APPROX 

SING  I  NT  EG 

0.5 

9881 

9985 

1.  0 

3544 

3401 

1.  5 

1269 

1251 

2.  0 

454 

476 

2.  5 

162 

167 

3.0 

58 

71 

3.5 

21 

24 

Table  17.  Values  of  J-integrnl  as  a  Function  of  Azimuthal  Scatter 
ing  Angle:  Z~  =  1.0,  T  =  15.0,  Exponential  Surface 


PH  ISC 


APPROX 

SIXG  I XT  EG 

POWER  EX  PAX 

9881 

3735 

7753 

3544 

1997 

45*59 

12(59 

1071 

1393 

■  454 

404 

551 

102 

200 

113 

58 

DC 

10 

21 

38 

2 

21 


35 


39 


2.  5. a  GAUSSIAN  RESULTS 


Examination  of  Tables  l  through  8  shows  that  the  SUM  representation  for  J 
(Gaussian  surface)  is  very  poor  for  E  »  l  and  becomes  more  and  more  accurate 
as  E  decreases.  For  E  »  1,  the  asymptotic  expression  APPROX  and  the  integral 
representation  for  J  agree  to  within  1  percent  for  all  PHISC.  As  E  decreases, 
the  asymptotic  representation  and  the  integral  representation  disagree  more  and 
more  with  each  other.  Finally  at  E  =  1,  the  asymptotic  representation  is  very 
inaccurate,  and  so  is  the  integral  representation  because  of  numerical  inaccuracies 
in  evaluating  INTEG.  The  problem  here  is  that  as  E  decreases,  the  decay  of  the 
integrand  as  r  increases  is  slow,  and  there  are  many  "c -dilations  o»  the  integrand 
(see  Appendix  B),  At  E  =  1,  and  E  <  1,  the  SUM  representation  of  J  is  the  only 
accurate  representation. 

2.  5.  3  EXPONENTIAL  RESULTS 

A  study  of  Tables  1 1  through  18  reveals  several  facts  about  the  behavior  of  J 
(and  hence  a9)  for  exponentially-distributed  surfaces.  For  large  Rayleigh  param¬ 
eters  E  »  1,  the  asymptotic  evaluation  of  the  Integral  for  J,  APPROX  and  the  inte¬ 
gral  representation  for  J,  SING  INTEG  agree  to  within  1  percent  for  all  azimuthal 
scattering  angles,  PHISC.  As  E  decreases,  APPROX  and  SING  INTEC  disagree 
more  and  more  with  each  other  until  at  E  ■  1.  the  asymptotic  representation  Is 
very  inaceu'-ate.  As  opposed  to  the  Gaussian  case,  for  the  exponentially -dis¬ 
tributed  surface,  the  integral  representation  for  J  is  accurate  for  all  values  of  the 
Rayleigh  parameter.  At  E  *  1,  the  asymptotic  representation  is  not  accurate, 
whereas  there  is  better  agreement  between  the  integral  representation  and  the 
power-series  expansion  representation.  For  E  *  0.  1,  the  asymptotic  representa¬ 
tion  is  extremely  inaccurate,  but  there  is  very  good  agreement  between  the  inte¬ 
gral  representation  and  the  power-?  Ties  expansion  representation. 

2.5.4  SUMMARY  OF  RESULTi 

To  summarize,  for  the  Gaussian-distributed  surface  heights,  the  asymptotic 
representation  and  the  integral  representation  for  J  are  accurate  for  large  Rayleigh 
parameters  E  »  1.  For  intermediate  Rayleigh  parameters,  the  sum  representation 
and  the  integral  representation  are  accurate.  For  small  Rayleigh  parameters 
E  <  1,  the  series  representation  is  the  only  representation  that  is  accurate.  On 
the  other  hand,  for  the  exponentially-distributed  surface  heights,  the  integral 
representation  for  J  is  accurate  for  all  Rayleigh  parameters.  The  asymptotic 
representation  is  accurate  for  large  Rayleigh  parameter  and  starts  to  deviate 
from  the  integral  representation  for  intermediate  Rayleigh  parameters.  For  small 
Rayleigh  parameters,  the  power-series  expansion  representation  is  accurate. 


3.  COMPARISON  OF  RESULTS  BETWEEN  STANDARD  GLISTENING 
SURFACE  AND  EXTENDED  GLISTENING  SURFACE 

0 

In  a  previous  report,  an  analysis  was  made  of  the  conditions  under  which  the 
conventional  definition  e.'  length  of  the  glistening  surface,  as  given  by  Beckmann 
and  Splzzichino1  (see  Appendix  C),  is  not  valid.  It  was  shown  that,  for  a  wide 
range  of  conditions,  significant  amounts  of  incoherent  scattered  power  can  be 
received  from  areas  beyond  the  conventional  length  of  the  glistening  surface.  The 
amount  of  additional  incoherent  power  not  included  under  the  conventional  efini- 
tion  depends  on  the  ratio  a! T,  except  for  the  case  where. both  transmitter  and 
receiver  are  very  close  to  the  surface.  It  was  shown  that  the  trends  in  the  be¬ 
havior  of  tne  glistening  surface  are  similar  for  vertical  or  horizontal  polarization, 
Gaussian-  or  exponential-surface  height  distributions,  and  for  different  signal 

•  g 

frequencies  in  the  S-band  to  L-band  region.  These  previous  investigations  exam¬ 
ined  the  effects  of  extending  the  length  of  the  glistening  surface,  while  the 
Beckmann  and  Spizzichino1  definition  of  glistening  surface  width  (VV)  was  used 
(see  Appendices  A  and  C).  The  diffuse  scattered  power  contribution  from  each 
increment  of  length  (St)  along  the  glistening  surface  is  obtained  by  multiplying  the 
product  of  the  centerline  value  of  a®  (#s  ~  0),  \V,  and  St  by  the  appropriate  azi¬ 
muthal  and  elevation  plane  antenna  power-pattern  distributions  (here,  -  azi¬ 
muthal  scattering  angle). 

In  this  report,  the  effects  of  extending  both  the  length  and  width  of  the  glisten¬ 
ing  surface  will  be  examined.  As  explained  in  Appendix  C,  the  length  is  extended 
by  integrating  the  product  of  a®  with  the  corresponding  elevation  and  azimuthal 
receiver  and  transmitter  gain  functions  over  the  entire  distance  between  the  two 
antennas.  The  width  of  the  glistening  surface  is  extended  by  integrating  the 
azimuthally  dependent  cross  section  o®  out  to  the  [Joint  where  the  azimuthal 
antenna  pattern  has  dropped  to  about  -45  dB  (average  sidelobe  level)  of  the  peak. 
Two  different  azimuthal  power  patterns  of  the  receiver  are  considered  in  these 
investigations:  a  monopulse  difference  pattern  (shown  in  Figure  2)  and  a  power 
pattern  that  has  no  null  at  boresight  ($  =  0.0).  This  second  power  pattern  is 

similar  to  the  original  monopulse  difference  pattern,  except  that  the  peak  at  about 
±3°  remains  constant  over  the  azimuthal  angle  range  -3°  s  6g  <3°. 

Figures  3  and  4  show  the  diffusely  scattered  power  (in  Watts)  as  a  function  of 

range  (in  nmi  and  km)  fo:  the  monopulse  receiving  antenna.  The  rough  surface  is 

uniform  with  a  correlation  length  T  =  100  m  and  three  variances  in  surface  height: 

2  2  2  2 

o  =  1  m  ,  10  m  ,  and  100  m  .  The  em  wavelength  A  =  0.275  m  and  the  complex 
dielectric  constant  of  the  surface  e  =  80  +  j9.0. 

Figure  3  is  a  triple  plot  of  diffuse  power  (Ppjpp)  vs  range  between  transmit¬ 
ter  and  receiver  where  the  Beckmann  and  Spizzichino*  definition  of  glistening 


Figure  2.  Azimuthal  Difference  and  Azimuthal  Sum  Pattern  and  Elevation  Pattern 
of  DABS  Antenna 


surface  length  and  width  was  used.  It  may  be  noted  from  this  figure  that  the  cut- 
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off  in  Pqjpp  at  short  ranges  for  o  -  1  m  and  a  =10  m  disappears  as  the  vari¬ 
ance  increasf  3  to  <P"  -  100  m^.  The  cutoffs  in  Ppjpp  at  short  ranges  is  due  to 
the  vanishing  of  the  glistening  surface  according  to  the  Beckmann  and  Spizzichino* 
definition.  In  Figure  4,  the  extended  length  and  width  definition  of  glistening  sur¬ 
face  was  used  to  calculate  Ppfpp  vs.  range.  It  may  be  noted  that  extending  the 
length  and  width  of  the  glistening  surface  removes  the  cutoffs  in  PQjpp  vs  range, 

so  that  the  glistening  surface  exists  at  all  ranges  for  all  values  of  the  variance 

2  2 

Also,  comparing  Figure  -i  with  Figure  4,  it  may  be  noted  that  for  a  *  10  m  and 
2  2 

o“  =  100  m“,  extending  the  length  and  width  of  the  glistening  surface  causes  an 

2  2 

increase  in  P.,,,.,.  for  all  ranges.  On  the  other  hand,  for  a  -  1  m  ,  P_IT_„  for 
DU'  I'  D1FF 

the  extended  definition  is  greater  than  Pjjjpp  for  the  standard  definition  only  for 
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Figure  3.  Diffuse  Power  Behavior  for  Azimuthal  Difference  Power., 
Pattern;  Standard  Glistening  -Surface,  T  -  100  m,  a-  :  i,  10,  100  m“ 
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Figure  4.  Diffuse  Power  Behavior  for  ^zimufhrl  Difference  Power 
Pattern:  Full  Integration,  T  =  100  in,  o~  ■  1,  10,  100  in- 


\ 


2  2 

short  ranges.  For  long  ranges,  ct  =  !  m  results  in  a  sharp  dropoff  of  a*  with 
increasing  azimuthal  angle  <j>s,  so  that  Ppjpp  for  the  extended  definition  Is  less 
than  Ppjpp  for  the  standard  definition  (which,  as  described  in  Appendix  C,  is 
obtained  by  multiplying  cr°  (bg  =  0)  by  the  classical  width  of  ihe  glistening  surface). 

All  the  parameters  in  Figures  5  and  6  are  the  same  as  in  Figures  3  and  4, 
except  for  the  receiving  antenna  pattern  that  is  now  centrally  uniform.  Figure  5 
corresponds  to  the  standard  definition  of  glistening  surface  and  Figure  6  corre¬ 
sponds  to  the  extended  length  and  width  definition.  By  making  comparisons  be¬ 
tween  these  four  figures,  it  may  be  noted  that  the  differences  in  Pfjjpp  between 
the  extended  definition  and  the  standard  definition  are,  in  general,  the  same  for 

the  modified  antenna  pattern  as  for  the  monopulse  pattern.  There  is  the  disappear- 
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ance  of  the  cutoffs  at  short  ranges  for  a  «  1  m  and  a  =  10  m  ,  and  there  is 
always  an  increase  in  PDjpp  when  the  glistening  surface  is  extended.  The  one 
notable  difference  is  that  for  the  modified  antenna  pattern,  extending  the  glisten- 
ing  surface  causes  Ppjpp  to  increase  at  long  ranges  for  a  *  1  m  .  For  the 
original  monopulse  pattern,  Ppjpp  decreased  as  the  glistening  surface  was 
extended.  This  difference  in  behavior  may  be  attributed  to  the  fact  that  removing 
the  null  at  boresight  in  the  receiving  pattern  causes  an  increase  in  PDIpp  (because 

a0  is  a  maximum  at  <$  .  =  0°),  which  compensates  for  the  falloff  in  a°  as  6 

s  s 

increases. 

To  conclude  this  section  we  will  compare  the  results  of  the  two  models  with 
data  taken  at  the  DABS  site  by  Lincoln  Laboratory  personnel.  This  data  consisted 
of  the  coherent  power  in  the  sum  channel  and  the  boresight  angular  uncertainty  in 
the  difference  channel  of  a  monopulse  receiver  located  on  the  ground.  The  trans¬ 
mitter  was  located  on  an  aircraft.  The  data  consisting  of  monopulse  boresight 
angular  uncertainty  vs  range  was  processed  by  removing  the  bias  in  boresight 
pointing  error  due  to  electrical  imbalances  in  the  system  and  then  converting  the 
azimuthal  angular  uncertainty  into  an  equivalent  diffuse  power  using  the  equations 
and  system  parameters  given  in  Appendix  A.  The  crosses  in  Figure  7  show  the 
results  of  the  processed  DABS  data. 

678  18  19 

In  Appendix  A,  and  in  previous  work  '  *  *  ’  it  was  discussed  how  statis¬ 

tical  estimation  theory  was  applied  to  digitized  terrain  maps  to  determine  param¬ 
eters  that  may  be  used  in  Vhe  em  scattering  models.  The  techniques  were  applied 
to  terrain  at  an  eastern  Massachusetts  site  (DABS).  A  data  base  of  topographic 
elevations  for  this  area  islnvailnbte  at  the  Electromagnetic  Compatibility  Analysis 
Outer  (L'CAC)  in  Annapolis,  Md.  This  was  prepared  from  digitized  terrain  maps 
supplied  bv  the  Defense  Mapping  Agency  (DMA),  The  area  of  interest  is  divided 
into  rectangular  cells,  each)  with  sides  about  2  km.  The  statistical  analysis  is 
then  ooplied  to  the  individual  cells. 
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Figure  5.  Diffuse  Power*  Behavior  for  a  Centrally-uniform  Power 
Pattern;  Full  Integration,  T  =  100  tn,  =  1,  10,  100  m2 
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Figure  7.  Comparison  Between  Standard  Glistening  Surface, 
Full  Integration,  and  DABS  Data:  Diffuse  Scattered  Power 
vs  Range 


The  results  of  the  terrain  data  base  analysis  are  recorded  on  a  computer  tape 
for  use  with  the  computer  program  for  the  electromagnetic  analysis  (see  Appen¬ 
dix  A).  The  dashed  line  in  Figure  7  shows  the  diffuse  power  predicted  by  the  com¬ 
puter  model  when  the  extended  length  and  width  glistening  surface  is  used.  The 
solid  line  shows  the  diffuse  power  predicted  by  the  computer  model  when  the 
Beckmann  and  Spizzichino1  definition  of  glistening  surface  is  used.  It  may  be 
noted  from  this  figure  that  the  extended  glistening  surface  model  shows  remark¬ 
able  agreement  with  experimental  data  for  all  ranges  between  transmitter  and 
receiver.  On  the  other  hand,  the  standard  glistening  surface  model  shows  poor 
agreement  with  the  experimental  data  particularly  at  short  ranges. 


4.  MOISTURE  EFFECTS 

In  this  section  we  are  concerned  with  the  behavior  of  the  scattering  surface 

when  there  is  a  range  of  moisture  levels  present.  The  moisture  is  represented 

by  assigning  different  dielectric  constants  to  the  surface.  This  effect  is  examined 

for  various  roughness  levels.  Both  definitions  of  the  glistening  surface  are  used 

and  their  results  compared  (Appendix  A  and  Appendix  C). 

7  19 

In  previous  studies  *  we  have  used  the  appropriate  dielectric  constants  for 
the  terrain  as  described  by  the  geologic  coding  of  the  DABS  site  data  base.  When 

O 

we  turned  to  the  question  of  glistening  surface  limitations,  we  used  both  the  DABS 
values  and  a  number  of  cases  where  the  entire  surface  was  assumed  to  be  of  uni¬ 
form  roughness  and  to  be  cultivated  terrain  with  complex  dielectric  constant, 
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=  80  +  j  9.0.  In  this  present  instance,  we  use  three  soil  descriptors  reflecting 
different  moisture  levels  and  consider  how  the  relative  roughness  affects  the  scat¬ 
tering.  For  dry  desert -type  sandy  soil  we  used  e  ,(1)  »  2  +  j  1.  62,  for  ordinary 

r  q 

dry  ground  et.(2)  =  4  +  j  0.006,  and  for  moist  soil  e  r<3 )  =  30  +  j  0.  6. 

Figure  8  shows  several  aspects  of  the  effect  of  moisture  content  on  the  scat¬ 
tering.  The  two  upper  illustrations  and  the  lower  left-hand  side  illustration  show 

the  progression  in  behavior  of  the  range-dependent  diffuse  scattered  power  when 
2  2 

a  *  10  m  and  T  is  varied  from  1,0  m  <  T<  500  m.  The  three  cases  are  for  the 
extended  definition  of  glistening  surface.  The  final  illustration  of  the  figure  shows 
the  solution  for  the  standard  definition  of  the  glistening  surface  with  T  =  500  m. 

If  we  examine  the  three  dielectric  constants  used  here  we  see  that,  although 
the  relative  change  in  the  imaginary  components  is  greater  than  that  for  the  real 
parts,  the  absolute  levels  are  smaller  in  magnitude  than  the  real  contributions 
and  thus  the  differences  in  the  diffuse  power  follow  the  real  part  of  the  dielectric 
constant  for  the  particular  set  of  soil  types  investigated  in  this  study. 
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Figure  8.  Soil  and  Moisture  Effects  on  Diffuse  Power:  Full  Integration  - 
o^  =  10  m^,  T  =  1,  100,  500  m;  Standard  Surface  -  =  10  m  ,  T  =  500  m 
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The  next  factor  to  note  is  that,  for  the  extended  surface  case,  as  the  value 
of  T  increases  from  T  *  1  m  to  T  =  500  m  for  the  extended  surface  case,  several 
differences  in  behavior  can  be  seen.  First,  the  respective  set;-  of  curves  show 
increased  rates  of  dropoff  with  increasing  range.  Next,  the  m  gnitudes  pass 
through  a  peak  at  the  intermediate  a/T  ratio.  Finally,  the  effect  of  moisture  con¬ 
tent  on  the  diffuse  power  increases  as  T  increases. 

The  standard  glistening-surface  results  for  T  =  500  m  differ  from  equivalent 

results  with  the  extended  glistening  surface  definition  as  is  expected.  For  the 

standard  case  there  is  the  typical  cutoff  at  28  nmi  and  the  overprediction  of  the 

8 

diffuse  power  at  long  ranges.  Here,  the  interest  is  in  the  effect  of  moisture 
level.  For  the  assumptions  of  the  standard  glistening-surface  model,  the  three 
moisture  curves  reveal  only  slight  variations  in  diffuse  power.  In  contrast,  for 
the  extended  surface  model  there  is  considerable  moisture  dependence  for  the 
scattered  power  (differences  of  as  much  as  an  order  of  magnitude).  This  is  a 
rather  significant  finding.  One  would  ordinarily  expect  that  the  relative  spread  in 
the  scattered  power  would  be  equivalent  in  the  two  models.  Some  analysis  of  the 
situation  will  be  useful  in  explaining  this  behavior. 

The  analysis  involves  how  the  moisture  levels  enter  the  calculation  and  hew 
the  two  versions  use  the  moisture  dependent  terms  to  arrive  at  the  diffuse  power 
levels  scattered  by  the  surface.  These  topics  are  addressed  in  the  discussions  of 
the  differences  between  the  two  scattering  models  contaircu  in  Appendix  C.  Thpre 
it  is  shown  that  the  dielectric  constant  is  introduced  as  a  term  in  the  Fresnel 
reflection  coefficients.  Then,  consideration  of  the  two  distinct  versions  of  the 
scattering  matrix  elements  indicates  how  the  difference  in  behavior  occurs.  In 
the  standard  model  the  scattering  elements  each  contain  only  one  of  the  Fresnel 
coefficients,  while  in  the  extended  version  each  element  depends  on  both  Fresnel 
coefficients.  In  addition,  for  the  standard  solution  the  final  result  tor  the  diffuse 
power  is  based  on  the  centerline  a’  variation  multiplied  by  the  glistening  surface 
width  at  each  point.  For  the  extended  solution  the  result  includes  the  scattering 
element  contribution  at  each  azimuthal  position.  This  introduces  a  wnole  range  of 
weighted,  moisture -dependent  terms  related  to  the  angular  variation. 

One  specific  aspect  of  the  particular  antenna  system  used  in  the  study  might 
suggest  that  these  results  are  only  valid  in  this  instance.  That  aspect  is  the 
receiving  antenna  azimuthal  power  pattern  (see  Figure  2),  The  pattern  contains  a 
centerline  null  that  removes  the  centerline  contribution  from  the  total  diffuse 
power  calculated  using  the  extended  surface  model.  The  question  is  whether  it  is 
the  aspect  that  allows  the  azimuthal  contributions  to  dominate  the  diffuse  power 
result  and,  consequently,  whether  there  would  be  moisture-related  differences 
between  the  two  models  if  alternative  azimuthal  power  patterns  were  considered. 
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To  answer  that  quest  on  the  original  power  pattern  was  replaced  by  one  having 
its  peak  at  the  centerline;  it  then  remains  at  ihat  level  over  a  considerable  frac¬ 
tion  of  the  entire  pattern.  When  the  diiiusc  power  was  calculated  for  the  new  power 

pattern  using  both  models,  differences  in  the  moisture  behavior  between  the  two 

2  2 

models  still  occurred.  Figure  9  shows  this  result  /or  the  case  where  o  =  10  m 
and  T  =  100  m.  For  the  two  antenna  azimuthal  power  patterns,  diffuse  power 
levels  are  presented  for  both  the  standard  surface  and  the  extended  one.  Close 
examination  of  these  results  indicates  that  the  second  power  pattern  results  in 
less  of  a  difference  in  moisture  effects  than  the  original.  Nevertheless,  even  with 
centerline  contributions,  the  extended  glistening-surface  solution  has  a  greater 
variation  with  moisture  level  than  does  the  standard  surface  solution. 


SOIL  AND  MOISTURE  CONTENT  EFFECTS  -  1 1 


»2«iom2  t*ioom 


MOV  (OIL 


EXTENDED  SURFACE  - MOIST  toil 


7  2 

Figure  9.  Soil  and  Moisture  Effects  on  Diffuse  Power,  a"  =  10  m  ,  T  =  100  m: 
Both  Antenna  Patters,  Standard  Surface,  and  Full  Integration  Models 
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Th"  final  area  of  the  nmisturc  analysis  is  a  summarization  of  the  offsets  for 
tiio  two  glisiening-surfnce  definitions  ns  ti  function  of  roughness.  This  compilation 
of  results  is  presented  in  Figure  10,  t  he  extended  surface  results  for  the  two 
moisture  extremes  are  shown  on  the  left  of  this  figure,  while  the  corresponding 
standard  glistening-surface  results  are  shown  on  the  right.  Each  graph  contains 
separate  curves  of  the  diffuse  power  behavior  at  three  antenna  sepa  ations  as  a 
function  of  a  roughness  factor,  ct  'T,  This  represents  a  highly -comprrssed  de¬ 
scription  of  the  results,  and  the  information  requires  detailed  discussion. 


SOIL  AND  MOISTUL’E  CONTENT  EFFECTS  III 

extended  surface  standard  surface 


(x/T  ct/T 

Figure  10.  Soil  and  Moisture  Effects  on  Diffuse  Power  as  a  Function  of  a/T 
Ratio  at  Three  Antenna  Separations:  Standard  Surface  and  Full  Integration 
Models  -  (5  nm i  - •;  30  nm i  — . - ;  50  nmi  . .  . . ) 


If  we  first  examine  the  two  surfaces  separately,  we  see  that  the  trends  for  the 
a/T  dependence  are. relatively  similar  for  the  three  separations.  The  behavior  is 
not  the  same,  though,  if  we  compare  the  results  for  the  two  definitions.  The  ex¬ 
tended  surface  results  in  larger  diffuse  power  levels  than  those  from  the  corre¬ 
sponding  standard  surface.  In  addition,  the  extended  surface  result  peaks  at  lower 


a/T  values  than  does  the  standard  case,  and  we  can  see  differences  in  behavior 
at  large  a/T  values.  The  two  definitions  produce  significantly  different  patterns 
in  this  range.  For  the  extended  case,  the  curves  at  the  three  antenna  separations 
all  tend  to  converge  for  large  a/T.  This  holds  for  both  moisture  levels.  On  the 
other  hand,  for  the  standard  case  the  two  sets  of  moisture  levels  result  in  paral¬ 
lel  behavior  for  large  a/T  values,  although  the  moist  surface  behavior  is  more 
complex.  For  the  standaru  case  an  additional  anomaly  is  apparent  at  the  opposite 
end  of  the  a/T  range.  There,  the  results  for  the  5-nmi  case  are  no  longer  sim- 
ila.  to  the  30-  and  50-nmi  cases.  This  is  a  result  of  the  vanishing  of  the  glisten¬ 
ing  surface  that  tends  to  occur  at  larger  ranges  as  the  a/T  ratio  decreases.  This 
subject  is  discussed  in  detail  in  our  previous  report,  but  the  trend  can  be  seen 

here  bv  consideration  of  the  standard  glistening  surface  cases  in  Figure  3.  This 
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figure  shows  the  change  in  glistening  surface  cutoff  for  a  =  10  m  as  T  changes 
from  T  =  100  m  to  T  =  500  m. 

Ther-»  is  no  simple  statement  that  generalizes  the  wide  range  of  moisture- 
related  effects  and  their  dependence  on  the  definition  of  glistening  surface  and 
scattering  model.  It  should  just  be  noted  that  there  is  a  strong  roughness  effect 
and  that  the  two  models  can  yield  quite  different  results  for  some  o/T  ranges. 

There  are  several  implications  for  radar  applications  in  these  results.  The 
first  point  is  that  knowledge  of  the  moisture  content  is  not  sufficient  for  absolute 
prediction  of  the  results  of  system  performance.  For  different  roughness  levels 
the  effect  of  moisture  can  be  quite  different.  The  second  point  is  that  moisture 
effects  based  on  standard  calculations  may  be  at  variance  with  the  predictions 
obtained  for  the  effects  of  moisture  on  diffuse  scattering  when  the  extended  surface 
is  used  in  the  analysis. 


5.  CONCLUSIONS 

In  this  report  we  have  examined  our  rough  surface  scattering  theory  in  terms 
of  its  underlying  principles  and  assumptions  and  established  that  the  assumed 
level  of  surface  roughness  has  a  highly  significant  effect  on  the  overall  solution. 
The  results  of  the  analysis  reflect  a  wide  range  of  topics.  In  this  section  we  will 
summarize  these  results  and  their  implications. 

The  first  area  is  thf  discussion  of  the  various  theoretical  formulations  and 
the  conditions  for  their  validity.  The  two  main  theories  for  scattering  are  the 
perturbation  theory,  which  requires  the  surface  irregularities  to  be  small  com¬ 
pared  to  a  wavelength  with  the  surface  slopes  less  than  unity,  and  the  physical 
optics  approximation  where  the  radii  of  curvature  of  the  irregularities  are  large 
compared  to  a  wavelength  and  the  surface  slopes  less  than  unity.  Additional 


factors  that  can  effect  the  results  are  multiple  scales  of  roughness  and  differences 
in  the  statistical  characterization  of  the  scattering  surface.  Among  the  more  sig¬ 
nificant  statistical  aspects  is  the  property  that  for  some  distributions,  uncorrelated 
surface  heights  are  independent  while  for  other  cases  this  condition  is  not  satis¬ 
fied.  In  the  present  discussions  conditions  for  physical  optics  are  always  assumed 
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to  apply.  Then  the  surface  cross  section  u°  =  J  depends  on  the  scattering 

matrix  and  the  statistical  form  J.  A  major  concern  of  this  study  is  to  examine 
the  agreement  and  limitations  for  the  alternative  J-represcntations  as  a  function 
of  Rayleigh  parameter,  £. 

For  Gaussian  surfaces  when  £  »  1,  the  summation  result  is  poor  while  the 
asymptotic  and  integral  results  agree.  In  this  case  the  agreement  is  independent 
.  of  azimuthal  scattering  angle.  These  solutions  diverge  as  £  decreases,  and  the 
summed  result  agrees  with  the  integral.  In  the  ranges  where  £  «  1,  the  summa¬ 
tion  result  has  become  the  only  usable  representation  because  the  behavior  of  the 
kernel  in  the  integral  formulation  makes  accurate  numerical  evaluation  of  the 
integral  intractable. 

For  exponential  surfaces  the  results  are  slightly  different.  In  that  instance, 
the  integral  form  is  accurate  over  the  entire  range  of  £  values.  The  asymptotic 
solution  agrees  with  the  integral  one  for  large  £,  with  the  summation  being  com¬ 
pletely  at  variance.  In  contrast,  the  £  <  1  the  asymptotic  solution  does  not  agree 
with  the  integral  result  while  the  power  series  expansion  solution  does.  Again, 
there  is  no  apparent  trend  related  to  the  magnitude  of  the  azimuthal  scattering 
angle. 

The  next  area  of  interest  is  that  of  the  concept  of  glistening  surface.  The 
extent  of  the  surface  and  the  inclusion  of  the  related  scattered  power  in  the 
theories  used  in  our  analyses  has  progressed  through  several  stages.  The 
differences  in  results  are  clearly  roughness  dependent. 

The  standard  glistening  surface  as  described  in  Beckmann  and  Spizzichino* 
generates  reasonable  results  for  a  wide  range  of  conditions,  particularly  for 
large  ct/T  ratios.  We  found  that  for  small  a/T  ratios,  though,  the  length  of  the 
glistening  surface  either  vanished  or  omitted  significant  portions  of  the  surface 
between  the  two  antennas.  Extension  of  the  glistening  surface  to  include  scatter¬ 
ing  contributions  from  the  entire  surface  between  the  antennas  has  shown  that  the 
power  from  these  excluded  areas  could  be  quite  significant  for  some  conditions. 
This  represented  the  first  modification  of  the  formalism. 

The  original  Beckmann  and  Spizzichino  solution  assumes  that  the  azimuthal 
variation  of  the  cross  section  is  negligible  and  determines  the  result  by  consider¬ 
ing  the  product  of  the  centerline  (4  =  0,  boresight  direction)  cr®  values  at  any 

point  along  the  length  of  the  glistening  surface  and  its  associated  width  at  that 
point.  In  the  initial  modification,  this  assumption  was  maintained.  The  newest 


modification  relaxes  this  constraint.  As  in  the  previous  case  the  integration  is 
along  the  entire  separation  between  antennas.  Here  the  concept  of  specific  surface 
width  is  removed,  and  the  integration  is  carried  out  over  the  azimuthal  extent  of 
the  power  pattern.  Also,  an  associated  azimuthally  varying  value  for  the  cross 
section  of  the  surface  is  included  in  the  theory.  It  was  found  that  for  small  a/T 
latios  the  standard  solution  approximations  result  in  over-prediction  at  large 
antenna  distances  and  underprediction  at  short  ranges.  For  larger  o/T  values  the 
results  tended  to  agree. 

For  the  above  results,  the  azimuthal  variation  ir.  a8  is  modified  by  the  corre¬ 
sponding  azimuthal  power  pattern  of  the  receiver.  This  suggested  considering  the 
effect  of  a  drastically  different  azimuthal  power  pattern.  Specifically,  instead  of 
a  monopulse  pattern,  a  pattern  was  chosen  that  has  an  essentially  constant  value 
over  the  typical  azimuthal  extent  of  significant  cross  section  values.  The  results 
followed  the  same  trends  as  before  except  that  the  full  integration  is  always 
greater  than  the  standard  glistening  surface  results.  In  the  case  Of  large  cr/T 
ratio,  the  centerline  null  in  the  original  pattern  resulted  in  relative  agreement  for 
the  two  models.  (The  centerline  null  in  the  pattern  balances  the  omitted  surface 
area.)  Here,  the  peak  centerline  value  adds  to  the  difference,  and  the  full  inte¬ 
gration  results  remain  higher.  The  final  point  with  respect  to  the  different  glis¬ 
tening  surface  results  is  that  the  full  integration  case  gives  excellent  agreement 
with  the  DABS  data  while  the  standard  surface  results  show  areas  of  poor  agree¬ 
ment  particularly  at  short  ranges. 

Once  we  established  the  effects  of  roughness  on  the  agreement  between 
theories,  the  final  topic  of  the  present  study  was  introduced.  This  is  the  effect  of 
roughness  on  the  results  when  the  moisture  content  of  the  surface  is  varied. 

A  number  of  results  were  found.  For  the  azimuthal  difference  power  pattern, 
a  decrease  in  the  value  of  a/T  resulted  in  increased  variation  between  the  succes¬ 
sive  moisture  level  values  for  scattered  power  as  a  function  of  antenna  separation. 
For  small  a/T  values  the  difference  between  the  moisture  level  results  of  the  full 
glistening  surface  and  the  standard  surface  is  considerable.  The  behavior  of  the 
diffuse  power  as  a  function  of  moisture  content  ard  surface  roughness  at  given 
antenna  separations  is  quite  dissimilar  for  the  two  definitions  of  glistening  surface. 
This  is  particularly  true  for  extremes  of  large  and  small  a/T. 

The  centerline  power  pattern  null  tends  to  emphasize  differences  between  re¬ 
sults  from  the  two  models,  but  the  introduction  of  the  alternative  azimuthal 
pattern  still  resulted  in  differ,  nces  in  the  effect  of  moisture. 

As  a  result  of  these  studies  of  roughness  dependent  effects  on  scattering  it  is 
clear  that  use  of  simplified  models  for  the  scattering  surface  can  lead  to  consid¬ 
erable  errors  in  the  prediction  of  the  scattered  power,  most  notably  for  small 
a/T  roughness  parameter  conditions.  In  addition,  the  effect  of  moisture  is  quite 


different  depending  on  the  glistening  surface  model  particularly  at  the  extremes 
of  either  large  or  small  o/T  ratio. 

Predictive  models  of  radar  performance  that  do  not  properly  account  for 
diffuse  scattered  power  can  lead  to  design  criteria  that  may  either  be  more 
expensive  then  necessary  or  result  in  poor  operational  performance.  The  mois¬ 
ture  dependence  on  roughness  also  suggests  that  there  can  be  inadequate  estima¬ 
tion  of  performance  variability,  unless  those  effects  are  properly  taken  into 
account. 
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Appendix  A 

Description  of  Rough  Surface  Scattering  Analysis 


In  this  appendix  vve  describe  the  details  of  the  surface  scattering  formalism 
in  terms  that  are  common  to  the  two  basic  versions  of  the  program  used  to  obtain 
the  results  found  in  this  report.  Additional  discussions  of  distinctions  between 
the  two  models  will  be  discussed  in  Appendix  C. 

The  common  elements  involve  the  calculation  of  the  specular  and  diffuse 
scattering  from  rough  terrain,  which  Incorporates  statistical  parameter  estima¬ 
tion  applied  to  digitized  terrain  data  bases.  The  terrain  is  divided  into  homo¬ 
geneous,  isotropic,  rectangular  cells,  that  are  each  characterized  by  a  mean 
height,  variance,  degree  of  correlation,  statistical  height  distribution,  and  an 
appropriate  dielectric  constant.  The  model  includes  spatial  inhomogeneities 
from  cell  to  cell,  multiple  specular  reflection  points,  and  global  and  local  shadow¬ 
ing  with  explicit  shadowing  functions  for  Gaussian  and  exponentially-distributed 
surface  heights.  The  antenna  power  patterns  of  both  the  transmitter  and  mono¬ 
pulse  receiver  are  included  in  the  model,  and  signal  processing  ltsses  can  be 
considered. 

The  programs  normally  use  system  parameters  associated  with  the  L-band 
Discrete  Address  Beacon  System  (DABS)  to  allow  comparisons  with  experiment. 
These  are  defined  in  Table  A1  and  Figure  8,  External  inputs  include  the  complex 
dielectric  constant  for  each  surface  area,  the  coordinates  of  the  monopulse 
receiver,  the  velocity  and  initial  and  final  positions  of  the  aircraft  containing  the 
transmitter,  and  a  parameter  to  control  the  effects  of  shadowing. 


Table  Al.  Experimental  Conditions  for  Discrete  Address  Beacon 
System  Tests 


Front-end  receiver  noise  figure 

3  dB 

Gain  of  monopulse  receiver  (sum  pattern) 
antenna 

22.5  dB 

Gain  of  transmitter  antenna 

4  dB 

Height  of  receiver  antenna 

101  m 

Height  of  transmitter  antenna 

1220  m 

Signal  polarization 

vertical 

Peak  transmitter  power 

350  W 

Pulse  length 

20  psec 

Azimuthal  beamwidth  (receiver) 

3® 

Wavelength 

0.275  m 

Transmission  line  loss  factor 

3  dB 

System  processi  'g  loss 

2  dB 

Normalized  pattern  slope 

1.5 

The  description  of  the  models  may  be  divided  into  two  major  topics.  These 
are  the  techniques  required  to  assign  appropriate  statistical  properties  to  the 
terrain  and  the  specifics  of  the  electromagnetic  formulation.  Each  of  these  as¬ 
pects  requires  some  discussion. 

There  are  several  surface  feature  contributions  in  the  models.  Analyses  of 
the  scattering  from  rough  surfaces  consider  the  surface  heights  in  the  region  as 
pairs  of  scattering  elements  and  in  most  cases,  including  this  study,  assume  that 
the  height  distribution  can  be  described  by  either  a  bivariate  Gaussian  or  exponen 
tial  probability  density.  These  two  bivariate  densities  have  the  forms: 

(Gaussian) 
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( [Exponential) 
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Both  forms  have  the  same  set  of  defining  parameters,  the  mean  height  the  var- 

•j 

irmee  a",  and  the  correlation  function  C. 

In  the  present  analysis,  either  the  form  of  the  density  function  and  its  assoc¬ 
iated  parameters  are  assigned  as  initial  conditions  of  the  problem  or,  if  a  specific 
site  is  being  analyzed,  the  probability  density  and  the  values  for  its  parameters 

are  determined  by  statistical  techniques.  The  procedures  for  describing  a  site 

A1  A2 

have  been  discussed  ir.  some  detail  in  earlier  work.  ’  A  brief  outline  will  be 
included  here.  Tile  theoretical  formulations  will  be  described,  followed  by  a  dis¬ 
cussion  of  their  application  to]  actual  terrain. 

The  starting  point  of  the  analysis  is  to  determine  appropriate  estimators  for 
the  parameters  of  a  Gaussian  or  exponential  density.  We  consider  that  we  have 
sample  height  values  from  the  population  of  the  region.  The  sample  mean  is  then 
used  as  the  estimator  for  the  population  mean  height  and  the  sample  variance  is 
used  as  the  estimator  for  the  population  variance.  A  more  complicated  develop¬ 
ment  is  required  to  establish  values  for  a  normalized  covariance  estimator.^2 
This  is  then  used  to  obtain  the  correlation  length,  T. 

The  next  step  is  to  decide:  which  density  is  more  appropriate.  Different  deci- 

J  ^3 

sion  models  have  been  developed.  One  such  approach  is  the  following.  If  the 
heights  in  a  region  are  described  by  a  particular  bivariate  form  when  considered 
in  pairs,  then  the  population  ajs  a  whole  should  be  described  by  a  higher-order 
multivariate  density  that  has  the  required  bivariate  marginal  density.  A  hypothesis 
testing  procedure  i.-;  then  applied  to  the  two  alternative  multivariate  densities  that 
describe  the  entire  collection  of  height  values  in  the  region  simultaneously.  This 
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simple  alternative  hypothesis  test  is  equivalent  '<'  1  r.  i •  - 1 •  •  > » : j probability 
criterion  where  it  is  equally  likely  that  either  density  is  the  appropriate  one. 

For  the  comparison  of  the  scattering  theory  with  the  experiment,  these  tech¬ 
niques  had  to  be  applied  to  the  terrain  at  the  eastern  Massachusetts  site.  A  data 
base  of  topographic  elevations  for  this  area  is  available  at  the  Fleet roniagnetic 
Compatibility  Analysis  Center  (FCAC)  in  Annapolis,  M  l.  This  was  prepared  from 
digitized  terrain  maps  supplied  by  the  Defense  Mapping  Agency  (DMA).  The  at  a 
of  interest  is  divided  into  rectangular  cells,  each  with  sides  of  about  2  km.  Kacli 
cell  is  further  subdivided  into  a  10  by  10  grid  of  points.  The  statistical  analysis 
is  then  applied  to  the  individual  cells. 

The  statistical  data  for  each  cell  has  been  recorded  on  a  computer  tape  for 
use  with  the  program  for  the  electromagnetic  analysis.  Fach  cell  is. represented 
by  seven  descriptors.  The  first  two  entries  are  the  (x,  y)  coordinates  for  the 
center  of  the  cell.  Next  is  the  geological  code  (dielectric  constant)  for  the  cell. 
(The  predominant  feature  is  woods;  there  are  a  number  of  cells  containing  clusters 
of  lakes  and  ponds  and  a  few  town  sites  with  associated  cleared  areas.)  This  is 
followed  by  the  mean  and  variance  of  the  heights  in  the  cell  and  the  estimated 
correlation  length,  T  (the  units  of  length  are  in  meters).  The  final  quantity  is  the 
result  of  the  hypothesis  test. 

The  trajectory  of  the  beacon  aircraft  is  incorporated  into  the  computer  pro¬ 
gram  and  at  each  range  point  for  which  a  calculation  is  to  be  made,  the  required 
cells  and  their  descriptors  are  then  identified.  These  results  or,  in  the  general 
case,  the  equivalent  set  of  input  parameters  are  then  used  in  the  electromagnetic 
analysis. 

The  calculation  of  the  electromagnetic  fields  has  two  distinct  elements.  First, 
the  total  coherent  electric  field  at  the  receiver  is  calculated  by  using  the 

sum  pattern  of  the  receiver  antenna  in  the  following  expression; 

KCOH  =  KT 


direct  path  electric  field  at  the  receiver, 

gain  of  receiver  in  direction  of  multipath  ray,  where  A  is  ttie 
angle  between  the  direct  ray  and  multipath  ray, 

attenuation  factor  affecting  coherent  reflected  wave  due  to  surface 
roughness, 


where 

E-i 


G.(A  )  = 


F. 

} 


1  + 


G.(«  )R. 
J  m  j 


ikAR 


F. 

3 


(A3) 


-  complex  Fresnel  plane  wave  reflection  coefficients, 


AR  =  difference  in  path  length  between  direct  ray  and  reflected  ray. 

The  summation  over  j  represents  all  possible  specular  reflections  (there  may  be 
more  than  one,  due  to  unevenness  in  terrain).  Here,  coherence  means  there  is  a 
known  phase  relationship  between  the  em  field  leaving  the  transmitter  and  that 
reaching  the  receiver. 

The  next  aspect  is  that  of  the  diffuse  power.  The  diffuse  power  in  the  mono* 
pulse  difference  channel  is  calculated  from  the  expression 

„  /PTXLOSSX  \  T r(  ° TR<<42)GR 

DIFF  ’  \  (4r)3  JJJ  \  R*R* 
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=  system  processing  losses 
=  transmitted  power, 

=  wavelength, 

=  gain  (power)  of  transmitter  in  azimuth  (isotropic  pattern), 
=  gain  of  receiver  in  azimuth  (difference  pattern,  Figure  8), 


=  gain  of  transmitter  i:.  elevation  (isotropic  pattern). 


=  gain  of  receiver  in  elevation  (Figure  8), 

=  elevation  angle  between  boresight  and  point  on  glistening  surface 
for  transmitter, 

=  elevation  angle  between  boresight  and  point  on  glistening  surface 
for  receiver, 

=  range  between  transmitter  and  point  on  glistening  surface, 

=  range  between  receiver  and  poin.  on  glistening  surface, 

=  element  of  area  of  glistening  surface  which  is  illuminated  by 
beacon, 

=  azimuthal  angle  between  boresight  and  point  on  glistening  surface 
for  receiver,  and 


=  azimuthal  angle  between  boresight  and  point  on  glistening  surface 
"  for  transmitter. 

The  diffuse  power  integral  contains  an  expression  for  the  normalized  average 

a  4 

bistatic  rough-surface  cross-section,  a®,  which  comes  from  Ruck  et  al.  The 
expressions  derived  by  Ruck  are  quite  general  and  highly  complicated.  The  dif¬ 
ferent  forms  of  a”  used  in  the  models  will  be  discussed  in  Appendix  C.  Here  we 
will  sint ply  comment  that 


■  |j) 


pq' 


J  s 


where  ^  represents  the  scattering  matrix  and  S  introduces  the  local  shadowing 
contribution.  The  term  J  is  related  to  the  surface  height  distributions  and  the 
surface  slopes.  For  a  Gaussian  surface 


and  for  an  exponential  surface  height  probability  density 


where 

f  =  sin  0.  -  sin  B  cos  <4  , 
x  1  s  s’ 

f  -  sin  <t>  sin  0  , 
y  vs  s’ 

{•  =  -cos  0.  -  cos  0  , 

Z  1  S’ 

<j>s  -  azimuthal  scattering  angle  , 

0^  =  angle  of  incidence  (with  respect  to  surface  normal)  ,  and 
0s  =  angle  of  scattering  (with  respect  to  surface  normal). 

The  final  aspects  of  the  model  is  the  azimuthal  angle  error.  To  calculate 
this,  we  assume  that  the  spectral  width  of  the  diffuse  multipatn  is  narrow  compared 
to  the  bandwidth  of  tile  receiver/processor,  and  that  both  noise  power  and  diffuse 
multipath  power  are  Rayleigh  distributed.  For  the  conditions  of  the  DABS  system, 

A4.  Ruck,  G.T.,  Barrick,  D.  E. ,  Stuart,  W.  D. ,  and  Krichbaum,  C.K.  (1970) 
Radar  Cross  Section  Handbook,  Vol.  2,  Plenum  Press,  New  York. 
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the  decorrelation  time  oi  the  diffuse  multipath  power  is  of  the  order  of  10  sec, 
and  the  interpulse  period  is  of  the  order  of  10  *’  sec.  Also,  for  tile  DABS  experi¬ 
ment  test  site,  the  spectral  width  of  the  diffuse  multipath  power  is  of  the  order  of 

4 

100  Hz,  and  the  bandwidth  of  the  receiver  processor  is  about  5  X  10  Hz.  This 
shows  why  even  narrowband  Doppler  filtering  cannot  reduce  the  diffuse  multipath 
power  in  the  radar  resolution  cell  containing  the  target.  Under  these  conditions, 
the  total  amount  of  noiselike  interference  Nj  in  the  radar  resolution  cell  containing 
the  signal  is  given  by 


P  + 

t  DIFF 


(A  5) 


where 

N  =  noise  power  from  environment  plus  receiver. 

Tlie  error,  Og,  in  azimuthal  angle  pointing  accuracy  is  given  by  the  expression 
of  Barton  and  Ward^ 


/  °B 

\  km  /2STIR 


(A  6) 


where 

flg  =  azimuthal  beamwidth, 

STIR  =  Pcoh/^I  =  siSnal  to  interference  ratio  in  the  difference  channel, 
Pcoh  =  coherent  power  in  sum  channel,  and 
km  =  normalized  pattern  slope. 


A“>.  Barton,  D.  K. ,  and  Ward,  II.  R,  ( 1 9 f>£>)  Handbook  of  Radar  Measurement, 
Prentice-Hall  Inc. ,  Knglewood  Cliffs,  New  Jersey. 
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Appendix  B 


Derivation  of  Expressions  for  Normalized  Cross  Sections 

of  Rough  Surfaces 


In  this  appendix,  the  expressions  for  the  normalized  cross  section  o'  of  a 
rough  surface  are  derived  for  bivariate  Gaussian-distributed  surface  heights  and 
also  for  bivariate  exponentially-distributed  surface  heights.  The  approximation 
of  physical  optics  is  assumed  to  apply.  Under  the  assumptions  of  physical  optics 
the  expression  for  o°  may  be  written 

*'  ’  lSpq|2jS  • 

The  scattering  matrix,  |ppql  and  the  shadowing  function,  S  are  discussed  in 
Appendix  A  and  rippendix  C.  Here  the  discussion  will  be  limited  to  the  final  term 
in  the  expression  for  a°,  the  J  factor,  in  Appendix  A  we  discuss,  for  both  dis¬ 
tributions,  the  particular  form  for  J  that  applies  for  large  values  of  the  Rayleigh 
roughness  parameter  [E  =  (27tct/\)(cos  0.  +  cos  9g)].  Here,  the  discussion  is 
broadened  to  other  regimes  and  the  form  for  the  J  term  is  derived  for  large 
Rayleigh  parameter,  Z  »  1;  intermediate  Rayleigh  parameter,  Z  *  1;  and  small 
Rayleigh  parameter,  Z  <  1.  It  should  be  mentioned  at  this  point  that  when  the 
Rayleigh  parameter  is  large,  a  formulation  of  the  rough  surface  scattering  prob¬ 
lem,  which  is  more  complete  than  the  physical  optics  approach,  shows  that 


1! ' 

multiple  scattering  bwomi'S  important In  particular,  DeSanto  and  Shislia 
showed,  by  numerically  solving  an  integral  equation,  that  tot  a  Gaussian  distrib¬ 
uted  surface  w  ith  large  surface  roughness  (large  Rayleigh  paiameter)  and  small 
correlation  length,  single  scattering  grossly  underestimated  the  strength  of  the 

average  scattered  field.  Recently,  for  a  rough  surface  having  a  correlation  length 

jj2 

ntueh  smaller  than  the  en  wavelength,  Brown  has  shown  that  an  analysis  may  bo 
used  based  on  a  subsi.iute  surface  which  gives  rise  to  the  same  average  scattered 
field  as  the  true  surface.  It  is  demonstrated  that  the  average  scattered  field  de¬ 
pends  upon  the  number  of  interacting  areas  and  the  surface  roughness.  There  is 
a  specific  number  of  interacting  areas  that  dominate  the  average  scattered  field 
for  a  given  range  of  Rayleigh  parameter  E.  The  number  of  interacting  areas 
increases  as  the  Ravleigh  parameter  increases. 

In  the  present  report,  no  multiple  scattering  is  considered  and  the  approxima¬ 
tions  of  physical  optics  are  used  to  derive  the  expressions  for  o6  when  the  surface 
heights  have  different  distribution  functions.  Different  expressions  are  obtained 
for  n0  as  a  function  of  the  value  of  the  Rayleigh  parameter.  On  Section  2  of  this 
report,  each  of  the  expressions  for  a’  is  evaluated  over  wide  range  of  the 
Rayh’igh  parameter  to  assess  tie-  regions  of  Validity  for  the  different  relations 
and  the  I'xtent  of  the  disagreement  among  the  o'  values. )  The  differences  in  <j- 
are  considered  to  he  related  to  the  differences  in  the  J -values  for  the  three  E- 
regimes.  The  generalized  integral  formulation  for  J  is  given  bv 


Ail 

f  dx2  /  d>  1  / 


<Jvf 


rxp 


i[vx(x2  -  V  +  vv(y2 


(Bl) 


where 

A  =  4XY 

vx  =  (ir){sin  °i  * si"  % oos  V 


Bl.  DeSanto,  J.A.,  and  Shisha,  O.  (1974)  Numerical  solution  of  a  singular  inte¬ 
gral  equation  in  r  andom  rough  surface  scattering  theorv,  J.  Comp.  Phvs. 
la<  \o.  21:28(1.  ~ 

B2.  Brow  n,  G.  S.  (1982)  New  results  on  coherent  scattering  from  randomly  rough 
conducting  surfaces,  IP KK  Trans.  Antennas  Propag.  (to  be  published). 


vy  =  ^(sin  9g  sin  ^g)  , 

0.  =  angle  of  incidence  with  respect  to  surface  normal  , 

0g  =  angle  of  scattering  with  respect  to  surface  normal  , 

=  azimuthal  scattering  angle  , 

X  *  marginal  univariate  characteristic  function  , 

Xg  *  bivariate  characteristic  function  , 

X  =  em  wavelength  , 

2X  =  length  of  rough  surface  ,  and 
2Y  =  width  of  rough  surface  . 

This  expression  for  J  involves  no  assumptions  about  the  surface  height  distribution 
or  the  size  of  the  Rayleigh  parameter  (the  only  restrictions  are  that  the  radius  of 
curvature  of  the  irregularities  on  the  surface  be  large  compared  to  a  wavelength 
and  that  the  surface  slopes  are  small,  cr/T  «  1). 

The  next  step  is  to  evaluate  the  integral  for  specific  cases.  We  first  note  that 
the  bivariate  characteristic  function  (the  Fourier  transform  of  the  bivariate  sur¬ 
face  height  distribution  function)  is  a  function  of  the  correlation  coefficient  c(t) 
where  r  is  the  separation  of  two  points  on  the  surface 

t  -  y/fxj  -  x2)2  +  (yr  -  y2)2  . 

2  2 

The  following  results  are  for  the  Gaussian  form,  c(t)  =  exp  [-f  /T  ).  The  corre¬ 
lation  coefficient  has  the  following  properties:  c(0)  *  1,  c(v  —  «)  —  0,  and 
0  <  |c |  <  1.  This  behavior  will  be  used  in  the  specific  deviations. 

The  first  case  that  we  will  consider  is  the  form  for  the  Gaussian  surface- 
height  distribution.  The  three  different  ranges  of  £  are  covered  in  succession, 

B3 

A  discussion  of  this  case  is  presented  in  the  book  by  Beckmann  and  Spizzichino. 

The  present  development  is  given  to  clarify  the  comparisons  in  Section  2  and  to 
show  the  similarities  and  differences  between  these  results  and  those  for  the 
exponential  distribution. 

Tor  the  bivariate  Gaussian  surface-height  distribution  function  as  given  in 
Appendix  A,  the  corresponding  characteristic  function  is  given  by 

X2  '  [-E2d  -  r(r)) 

B3.  Berkmann,  P.  ,  and  Spizzichino,  A.  (lt*n:j)  the  Scattering  of  I T’ct rcmngnet ic 
Waves  I-'rnm  Rough  Surfaces,  Macmillan  Co.  .  Yew  York. 


and  (ho  corresponding  univariate  characteristic  function  is  given  by 
X  =  exp  [  ”/2]  . 

For  the  Gaussian  distributed  surface  heights 
\9(t  -  «)  -  x  X* 

that  is,  decorrelation  leads  to  statistical  independence.  Because  of  this  fact  the 
previous  integral  for  J  may  be  readily  reduced  to  the  following  form: 


2 

J  =  f  do(vxvT)[X2  ‘  X  X<!)  "  JT  1 


(B2) 


hherevxy  =  vj  +  vj. 

If  the  Rayleigh  parameter  E  »  1,  a  steepest  descent  evaluation  of  the  integral 
shows  that  most  of  the  contribution  comes  from  r  ~  0.  Then 
c(t)  =  exp ( -f  “  T“|  =  (1  -  T“/T“i  =>  '  exp  (-E2 t“/T”)  and  x  X*  *  0  so  that 


J  =  Trf  J0{\yT)  eXp  [*-2t2't2|  ** 


and  after  integration. 
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The  above  expression  is  accurate  for  large  E.  On  the  other  hand,  the  most 
accurate  representation  for  J  for  intermediate  values  of  E  is  given  by  tile  integral 
in  Kq.  (B2).  When  the  Rayleigh  parameter  is  small  (E  <  1),  the  following  series 
representation  should  be  used 
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tvz  |  V  1  -  C  W  +  (c  -  1)Z2  j 


-oo  -oo 


(exp  [-  c^w2  +  Z2)  '  ]}. 


Next,  we  cnange  to  cylindrical  coordinates,  where  w  =  r  sin  9  and  Zj  *  r  cos  9 
so  that 


/ - T  °r  /■*  -iv_r[  vd-c2  sin9+(c-l)cosS)  -c,r 

X2  ~  c  *  l  ~  c  J  r  dr  I  d6  e  e  1 


0  0 


Performing  the  integration  over  r  yields: 


*2  -  c" V  1  -  c2  y  d 0  -  +  iv^  (v 1  -  c2  sin  9  +  (c  -  1)  cos 


\]“2 
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This  integral  can  be  evaluated  by  the  calculus  of  residues  if  we  introduce  the 
change  of  variables  Z  =  e’^,  so  the  Xj  becomes 


_  -4  iFvl  -  c 

'2  '  9 

*  Ai 


[Z2  +  (B/A)Z  +  (D/A)]2  . 


=  vz  -i(I  -  c)  +  vd-c2  , 


B  =  -2c  ^  ,  and 


r  i 

1>  -  V  [i(l  - 


r)  t  VI  -  c“ 


and  lie-  contour  1st  a  round  the  unit  circle  in  the  complex  /.-plane.  Conchy's 
integral  formula  men  leads  to  the  result 


'  V/ 


X2  =  (1  +  (?/3)E“[l  -  c(t)H 


An  alternative  derivation  of  this  result  was  reported  by  Brown. 

The  marginal  univariate  density  function  corresponding  to  the  bivariate 
exponential  distribution  is  given  by 


p(Z9)  •  2c  v  1  -  c~ 


J  du  exp  -Cj  (u2  +  zl  j1 


27  VI  -  r 
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2rX^|z2|K1(c1|z2|) 


Thru,  thr  univa riate  characteristic  function  is  given  by 


r  ivz; 

=  J  p(Z2)  e 


and  hence, 


=  [l  +  (l/3)£“  ] 


Here,  it  may  br  noted  that  X2^T  -*  «)  *  \  \*  so  that  decorrelation  docs  not  imply 
statistical  independence. 

As  in  the  Gaussian  case,  if  the  Hayleigh  parameter  is  sufficiently  large 
(E  »  1),  then  \2  -  x\*  ~  \.,  and  the  contributions  to  the  four-fold  integral  [Eq.  (1)] 
for  J  come  primarily  from  the  region  where  T  s  0,  so  that 
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B4.  Br  own,  G.  S.  (1082)  S* *;i 1 1 # * r*i nu  from  n  rlnss  of  randomly  rough  surfaces, 
Radio  Sri,  (ro  bi*  published). 


/ 


which  reduces  to 


r  dr  . 


This  is  the  expression  given  by  Barrick  in  Buck  et  al.  ' 

On  the  other  hand,  if  the  Rayleigh  parameter  is  not  large  (£  <  1,  intermed¬ 
iate  Rayleigh  parameters),  then 


I  J 

J=  (tj)(t)  /dxi  /  dx2  /dy i  /dy2Q(x i.x2yi.y2)[x2 -xx*I  . 


-X  -X  -Y  -Y 


where 


Q<xx.  x2*  yr  y2> =  pxp  t‘<vx[xi  “  x2^  +  Vyi  '  y2* 


Brown  has  rewritten  this  expression  as  the  sum  of  two  integrals,  by  subtracting 
and  adding  the  term 

-3 /2 

X2(r  -  «,)  =  [1  +  (2/3)£2] 

which  follows  from  c(r  —  oo)  -  0.  Then,  J  =  Jp  +  Jg,  where 


X  X 


dD  J  dxl  /  dx2  /  dyl  /  dy2  Q<5Ci'X2.y1.y2){  t2  -  X2<t  -  ®>J 


-X  -X  -Y 


and 

Jc  =  /-~)(4XY)sinc2  (v  X)  sine2  (v  Y) 

a  x  y 

Here:  sinc(u)  =  (—  l~-  — ) .  Since  most  of  the  contribution  to  comes  from  the 
region  where  T  *  0,  one  may  write 

B5,  Ruck,  G.  T.,  Barrick,  D.  K. ,  Stuart,  W.  D. ,  and  Krichbaum,  C.  K.  (1970) 
Radar  Cross  Section  Handbook,  Vol.  2,  Plenum  Press,  Xevi  York. 
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Appendix  C 

Differentiating  Features  of  the  Two  Models 


As  was  discussed  in  Appendix  A,  many  of  the  details  for  the  two  models  are 
identical.  The  point  at  which  the  models  differ  is  in  the  calculation  of  the  diffuse 
scattered  power.  In  the  discussion  we  will  first  present  the  specifics  of  the  gen¬ 
eral  model  and  then  show  how  some  simplifying  assumptions  lead  to  an  alternative 

less  complex  formalism  that  includes  the  length  and  width  definitions  of  Beckmann 
Ci  C2 

and  Spizzichino.  Ruck  et  al  give  expressions  for  the  average  bistatic  rough- 
surface  cross-section  oq  under  the  following  four  assumptions:  (1)  the  radius  of 
curvature  of  the  surface  irregularities  is  larger  than  a  wavelength;  (2)  the  rough¬ 
ness  is  isotropic  in  both  surface  dimensions;  (3)  the  correlation  length  is  smaller 
than  either  the  X  or  Y  dimension  of  the  sample  subregion;  and  (4)  rr. ultiple  scat¬ 
tering  is  neglected.  The  expression  for  oq  has  been  presented  in  Appendix  A. 

Here  we  wish  to  discuss  the  scattering  matrix  element  that  is  used  in  the  deter¬ 
mination  of  the  cross  section  of  the  rough  surface. 

The  matrix  elements  for  linear  polarization  states  are 


VV 
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Note  that  er  is  the  relative  complex  dielectric  constant  of  the  surface,  the  sub¬ 
script  11  refers  to  the  10-field  in  the  plane  of  incidence,  and  the  Subscript  1  refers 
to  the  K-field  normal  to  the  plane  of  incidence.  The  remaining  angle  terms  are 
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Then,  neglecting  shadowing  effects  we  have  the  general  expression  for  the  surface 
cross  section 


*s>=  l%|2j  • 

The  diffuse  scattered  power  is  determined  by  integrating  the  product  of  this 
expression  and  the  corresponding  elevation  and  azimuthal  receiver  and  transmitter 
gain  factors  over  the  distance  between  the  two  antennas.  The  azimuthal  integration 
is  bounded  by  the  extent  of  the  azimuthal  difference  pattern.  Figure  10.  (For  a 

discussion  of  the  technique  used  to  introduce  these  bounds  into  the  formulation, 

C3 

see  Papa,  Lennon,  and  Taylor.  ) 

In  the  alternative  model  the  assumption  is  made  that  the  receiver  is  sufficiently 
far  from  the  transmitter  so  that  the  portion  of  the  "glistening  surface"  that  con¬ 
tributes  to  the  diffuse  multipath  is  a  long,  narrow  strip  extending  between  the 
transmitter  and  receiver.  This  assumption  allows  us  to  make  the  approximation 
that  the  azimuthal  scattering  angle  <£g  »  0°.  This  assumption  leads  to  considerable 
change  in  the  resultant  a°  calculation. 

In  the  scattering  matrix,  the  first  consideration  is  that  the  two  cross -polarized 
terms  are  now  zero.  Only  the  copolarized  elements  contribute  to  the  scattering. 

In  those  cases,  manipulation  of  the  relations  and  introducing  the  angle 

i  0.  +  0  \ 

a  -  ( — ^2 — — )  leads  to 

(1  +  cos  2o)R  (0.) 

0VV  =  (cos  0, +  cosO  )~  (vertical  polarization) 

1  S 

and 

(1  +  cos  2a)  R  (0.) 

PHH  "  (cos  0.-+  cos  0ST~  (horizontal  polarization)  . 


In  these  two  terms  it  should  be  noted  that  each  depends  on  only  one  Fresnel  reflec¬ 
tion  component  while  in  the  general  solution  each  includes  contributions  from  the 
two  Fresnel  coefficients.  Another  effect  of  this  assumption  is  seen  in  the  terms 

?  ,  ?  ,  and  S  .  These  reduce  to  £  =  sin  0.  -  sin  0  ,  ?  =0,  and 

x’  -y*  z  x  i  s’  y  ’ 

r  s  -cos  fl  -  cos  0  .  As  a  result  we  then  have 
z  is 
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for  a  Ga  ;ian  bivariate  surface-height  probability  density  function  (PDF)  and 
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Since  this  has  effectively  removed  any  azimuthal  variation  in  or®,  the  calcula¬ 
tion  of  the  diffuse  power  scattered  by  the  surface  and  received  at  the  antenna  is 
approached  differently.  In  this  model,  the  surface  of  integration  is  determined  by 
the  Beckmann  and  Spizzichino  definitions  of  glistening  surface  length  and  width. 

The  length  is  obtained  by  calculating  the  location  of  the  two  end  points  of  the  glis¬ 
tening  surface.  For  a  homogeneous  surface,  the  distances  from  the  end  points  to 

the  transmitter  and  receiver  (L^,  L^)  are  based  on  the  respective  heights  (H^.,  H^) 
and  a  roughness  criterion  (a/T),  where  is  the  surface  height  variance  and  T  is 
the  correlation  length.  Then, 

Lj  =  Ht  cot  <20o) 

L2  =  Ha  cot  (2pQ)  (Cl) 

where, 


tan  [io  =  2ct/T  . 


Then,  along  the  extent  of  the  surface  the  local  width  W  is  given  by 
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Here, 

D  =  total  ground  distance  from  transmitter  to  receiver, 

Xj  =  distance  from  transmitter  to  point  on  glistening  surface,  and 
X2  =  distance  from  receiver  to  point  on  glistening  surface. 

f!8 


Once  these  quantities  have  been  determined,  the  diffuse  scattered  power  contribu¬ 
tion  from  each  increment  of  length  along  the  glistening  surface  is  obtained  by 
multiplying  the  product  of  the  centerline  value  of  a0  and  W  by  the  appropriate 
azimuthal  and  elevation  plane  antenna  power-pattern  distributions. 

The  differences  in  the  various  results  described  in  Section  3  and  Section  4  of 

this  report  reflect  the  effects  of  using  the  two  alternative  models  in  our  analyses. 

C4 

One  additional  point  that  can  be  made  is  that  in  an  earlier  report  we  described 
an  attempt  to  go  beyond  the  Beckmann  and  Spizzichino  model  limits  without  resort¬ 
ing  to  the  complete  solution  described  here.  In  that  instance,  the  restriction  on 
the  surface  of  integration  based  on  the  length  of  the  glistening  surface  was  re¬ 
moved  and  the  calculation  was  made  over  the  entire  distance  between  antennas 
whether  or  not  that  exceeded  the  traditional  glistening  surface  length.  However, 

we  still  assumed  <4=0  and  the  standard  width  values  were  used, 
s 
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